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Abstract 

The gyrokinetic description of particle dynamics faces a basic difficulty when a special type of 
canonical variables is sought, i.e., the so-called gyrokinetic canonical variables. These are defined 
in such a way that two of them are respectively identified with the gyrophase-angle, describing the 
fast particle gyration motion around magnetic field lines, and its canonically conjugate momentum. 
In this paper we intend to discuss the conditions of existence for these variables. 

PACS numbers: 52.30.Gz,52.30.-q 
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I. INTRODUCTION 



The gyrokinetic description of particle dynamics concerns the representation - in terms 
of generally non-canonical variables - of the state of a single classical charged point-particle 
immersed in an electromagnetic (EM) field. In classical electrodynamics - for arbitrary EM 
fields and provided the EM self-force is neglected (according to the customary interpretation; 
see, however, the related discussion in Ref. - it is well known that this defines an 
Hamiltonian dynamical system. As such, locally in phase-space, it can always be represented 
in canonical form. Nevertheless, although well-known in the literature, the gyrokinetic 
problem presents a basic difficulty related specifically to the construction of a special set of 
canonical variables. As described below, this arises in connection with the construction of 
the so-called canonical gyrokinetic variables. In this paper we intend to show that, unless 
special symmetry conditions hold (in particular for the magnetic field), such variables do 
not exist. 

For a proper formulation of the problem it is important to stress that the gyrokinetic 
description of particle dynamics concerns in principle two possible viewpoints : 

A) an exact representation (of the particle state) realized by suitably prescribing an ap- 
propriate phase-space diffeomorphism [see below Eq.(J\)]; 

B) an approximate representation, obtained by means of a suitable asymptotic expansion. 

In the exact gyrokinetic treatment (Approach A) it is assumed that there exists a phase- 
space mapping (i.e., a diffeomorphism, which is assumed to exist at least in a suitable subset 
of phase-space) of the form 



between the Newtonian particle state (r, v) (with r and v denoting respectively the particle 
position and velocity) and the (real) gyrokinetic state (y',0'). The latter is defined in such 
a way that equations of motion for the transformed variables take the form 



By construction Y(y , t) and -F(y', t), to be identified with suitably smooth real functions of 
the gyrokinetic state, are both assumed independent of the variable 0'. Hence, its canonically 
conjugate momentum p^ is necessarily a first integral of motion. Here <$' and y' are the 




(1) 




(2) 
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so-called gyrokinetic (or guiding-center) variables, representing respectively an angle (the so- 
called gyrophase) which describes the particle gyration motion around the magnetic flux lines 
and an arbitrary 5-dimensional real vector, representing a reduced non- canonical gyrokinetic 
state. As an example, in particular, the vector y' may be identified with the non-canonical 
variables 

y = (r',f',P0') , (3) 

where r' and £' denote respectively the guiding-center position vector and an additional 
(independent) velocity-space gyrokinetic variable. 

However, unless the electromagnetic field is specially prescribed (i.e., for example, it is a 
constant), the exact gyrokinetic transformation (pT|) cannot generally be achieved. Neverthe- 
less, it is still possible - under suitable assumptions - to determine it in an approximate sense 
by means of an appropriate asymptotic approximation for the (canonical or non-canonical) 
particle state (Approach B). This is obtained by introducing an asymptotic expansion (i.e., 
a truncated perturbative expansion in terms of an appropriate infinitesimal dimensionless 
parameter e) of the form 

N 

(r, v) - (y', 0') = £ e k (y' k , ft) + o(e N+1 ) (4) 

k=0 

where the integer N > 1 (to be suitably prescribed) denotes the "order" of the asymptotic 
approximation. In particular, by assuming that the magnetic field in which the particle is 
immersed is suitably "intense" , the infinitesimal dimensionless parameter can be defined as 

e = r L /L « 1. (5) 

Here the notation is standard. Thus, L and n are respectively a characteristic scale length 
of the EM fields (to be suitably defined, see below) and the velocity-dependent particle 
Larmor radius = f^. In particular, all primed quantities are evaluated at the guiding- 
center position r', which requires that the diffeomorphism 

JV 

r ^ r' - J>V fc + (^ +1 ) (6) 

is assumed to exist. For example Q' = Q'(r',t), where Q' = ^ is the Larmor frequency 
and q, m, B are respectively the charge and mass of a point particle and the magnitude of 
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the magnetic field. Moreover, w' is the orthogonal component of the particle velocity to be 
evaluated - in a suitable reference frame - at the same position r'. 

The first author who systematically investigated the gyrokinetic problem, based on the 
explicit construction of an asymptotic expansion of the form (J3J) and was Alfven |2J who 
pointed out the existence of an adiabatic invariant, the magnetic moment // = ^ p'^, in the 
sense: 

|]n„'~0( e ). (7) 

After subsequent work which dealt with direct construction methods of gyrokinetic variables 
3, Q, E|, a significant step forward was made by Kruskal [6] who, first, established the 
consistency of the Alfven approach by proving, under suitable assumptions on the EM 
fields, that the magnetic moment can be constructed correct at any order in e in such a 
way that, denoting M' such a dynamical variable, it is an adiabatic invariant of order N, 
namely in the sense 

-^-lnM' ~ 0(e N ) (8) 
W = // + + ... + £ Vat- (9) 

A modern picture of the Hamiltonian formulation, which makes easier the formulation of 
ligher order perturbative theories, was given only later in terms of Lie-transform methods 
7, sl, 1^]]. However, it was only with the adoption of non-canonical Lie-transform methods 
9] that the approach was given a general formulation. As a motivation to his non-canonical 

on 

approach, Littlejohn |9J, H0[ pointed out what in his views was a critical point of purely 
canonical formulations such as previously developed Lie transform approaches s| , namely 
the ambiguity in the separation of the unperturbed and perturbed contributions in the 
Hamiltonian due the presence of the vector potential A in the canonical momenta. He 
showed that this difficulty can be circumvented by making use of suitable non-canonical 
variables which include the canonical pair (<f)' , p'^) . 

The possibility of constructing canonical gyrokinetic variables has relied, since, on two 
possible methods 

• the Darboux reduction algorithm, based first on the construction of a set of non- 



canonical gyrokinetic variables 111 ]; 
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• the direct construction of canonical gyrokinetic variables, either in terms of mixed- 
variables generating functions canonical Lie-transform methods ^_8, 12] or based 
on the adoption of the hybrid Hamilton variational principle ll[ 13, 14 ]. 

The first approach, and probably the most popular in the literature HQ, HQ, is the 
based on the use of Darboux theorem which allows, in principle, the construction of canoni- 
cal variables for an arbitrary differential 1-form. The canonical 1-form expressed in terms of 
the canonical variables is then obtained by applying recursively the so-called "Darboux re- 
duction algorithm" as pointed out by Littlejohn, which is obtained by a suitable combination 
of dynamical gauge and coordinate transformations. Nevertheless this approach leads to 
potential complications and ambiguities due to the fact that the so-called canonical gyroki- 
netic coordinates (see below) are field-related. Therefore it would be highly desirable to be 
able to construct gyrokinetic canonical variables which result independent of the magnetic 
field geometry. In a previous work Tessarotto and Nicolini, 2006 [19]) a possible solution to 
this problem has been pointed out by adopting superabundant canonical variables. Purpose 
of this work is - instead - to address the problem of the construction of essential gyrokinetic 
canonical variables. 



II. LAGRANGIAN APPROACH AND CANONICAL GYROKINETIC REDUC- 
TION 

Starting point for the application of the Darboux reduction method is the standard 
Lagrangian formulation for gyrokinetic particle dynamics, expressed in non-canonical gy- 
rokinetic variables. For definiteness, let us briefly recall its formulation. Let us assume for 
this purpose that the EM potentials $, A are analytic functions of e and hence can be 
represented in power series of e 

N 

$ = Y,e%(r,t), (10) 

i=-l 
N 

A = J2^Mr,t). (11) 

i=-l 

Here e is the infinitesimal dimensionless parameter defined above ([1]). In particular, the 
characteristic scale length L entering its definition is identified with the minimum of 
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the gradient-scale lengths for the perturbations of the EM potentials ($«, A*), namely as 
L < min (r^ ^§jr ^ , m i n f° r * = —1,N. In addition, denoting where 

b(r, £)= B(r, t)/B(r, t), the magnitudes of the particle velocity |v| and of the electric drift 
velocity vg = cE x b/5 are assumed of the same order, in the sense 

\v\ / K\ ~ o(l) (12) 

and consequently the parallel electric field is similarly ordered as 

b-E~o(l) (13) 

(condition of small parallel electric field). In validity of these hypotheses the construction of 
the standard gyrokinetic variables is well known and has been achieved by several authors 
(see for example [9]). For definiteness, let us identify the reduced gyrokinetic state y' with 
y' = (r', u', p'A . Here vl denotes the parallel velocity 



u 



> = b>. (V-V E ) , (14) 



where v' is the guiding-center velocity and V E = — J?b' x V'0' is l ne E x B-drift velocity 
evaluated at the guiding-center position. In this case - and in the presence of slowly varying 
EM fields - the fundamental Lagrangian differential 1-form expressed in terms of gyrokinetic 
variables reads 

dV = dt£'(y', v', r', 0', t) = dG' - d<f)'p'^ - dtH' (15) 
where dG' and H' are respectively the exchange term 

dG' = a(y',t) -dr', (16) 
a (y',t) = ^AV.*). ( 17 ) 



and the gyrokinetic Hamiltonian 



y v' 2 + + V(y',t). (18) 

Moreover, for definiteness, let us identify the reduced state y' with a suitable "effective" 
vector potential A*(y', t), i.e., it reads . Here, both the exchange term and gyrokinetic 
Hamiltonian, in particular the effective EM potentials ($*, A*) are expressed as functions 
only of the non-canonical reduced gyrokinetic state y' defined above ([3]). In the following the 



gyrokinetic differential 1-form dV will be considered either exactly prescribed (Approach A) 
or determined in terms of an asymptotic approximation of order o(e N+1 ), namely neglecting 
corrections of order o(e N+1 ) to dV (Approach B). 

Let us now seek a diffeomorphism, to be assumed at least locally defined in the relevant 
phase-space, of the form 

(r',M',p0/,0') -> [q a ,q' 2 ,p'i,P2,P^Aj , (19) 

where y' 1 , y' 2 , p x and p' 2 are assumed smooth real functions only of y' (and hence by definition 
as gyrokinetic variables). Provided the differential 1-form dG' when expressed in terms of 
[q >l , q' 2 , p' x , p' 2 , pp) takes the canonical form 

dG' = p' 1 dq n +p' 2 dq' 2 , (20) 

(canonical gyrokinetic reduction) the variables z'— (q' 1 , q' 2 , p'i, p'2, Pft , 0') are manifestly 
canonical gyrokinetic variables. In fact, it is immediate to prove that the Euler-Lagrange 
equations corresponding to the variational differential form ffT5]) expressed in the variables 
z' are canonical. 

Particle dynamics expressed in terms of the canonical variables z' denotes the so-called 
canonical gyrokinetic treatment (CGKT). The explicit construction of these variables has 



been first pointed out by Littlejohn [11] . adopting the so-called Darboux reduction method, 
by considering the vector potential A, and hence the associated magnetic field B (equilibrium 
magnetic field), as stationary. However, the proof - achieved in this way - of the local 
existence of the diffeomorphism ffT9]) and hence of the gyrokinetic canonical variables defined 
above z', is not generally applicable to general situations. In fact, to reach it in Ref. ll| 
it was assumed that the magnetic field admits, at least locally (in configuration space), a 
family of nested toroidal magnetic surfaces. 

This raises, therefore, the issue of the general validity of such a conclusion. In fact, 
the question is whether it applies only in the case of equilibrium magnetic fields which 
are symmetric, i.e., which possess at least one ignorable coordinate, or - at most - exhibit 
suitably small deviations from a symmetric equilibrium. In fact, it is well known that the 
proof of existence of smooth MHD equilibria with good magnetic surfaces (namely which 
admit a family of locally nested toroidal magnetic surfaces in a finite subset of configuration 
space) can only be achieved for symmetric equilibria (1) or at most for magnetic fields which 



are asymptotically close, in some sense, to equilibria of this type (2). As an example, in Ref. 
[3] to obtain the canonical variables with the Darboux reduction method, consistent with 
the requirement (2), it was assumed a magnetic fields almost axi-symmetric, i.e., allowing 
actually only infinitesimally small deviations from axi-symmetric toroidal geometry. 



III. ON THE EXISTENCE OF CANONICAL GYROKINETIC VARIABLES 

For definiteness let us pose, in this Section, the problem of the existence of the canonical 
gyrokinetic variables in the framework of the exact gyrokinetic formulation (Approach A). 
In order to solve the related problem let us analyze the conditions of validity of the canonical 
gyrokinetic reduction (1201) in the particular case in which there results a = a(y ) in dG' [see 
Eq. ffToT) ]. For this purpose let us seek a diffeomorphism 

r'^q'^'W^O/VV 3 ) (21) 

denoting in principle arbitrary real and gyrokinetic variables. These can be defined, in 
particular, in such a way that 

(rV,jy,0') -> (gVV 3 ,^,^) , (22) 

with £' to be suitably defined, is a phase-space diffeomorphism. Hence it follows that the 
differential 1-form dG' has necessarily the general form 

dG' = f'idq' 1 = dG" (23) 

where // = f[ (q', £',£><//) • The analysis of the conditions of validity of the dynamical reduc- 
tion - under which the differential 1-form dG" , as given by Eq. (|23[) . can be brought to its 
canonical form (1201) - is straightforward. Let us first establish the following lemma 
Lemma - Reduced form for dG" 

Let us assume that the real functions f- = f[ (q', (for i = 1,2,3): 

1) are suitably smooth (i.e., at least C^ 2 '); 

2) the set of gyrokinetic variables (q n , q' 2 , fx, f2,P<f>', (ft') are defined so that they all in- 
dependent; 

3) are defined so that for at least an index i (for i = 1,2,3) there results 
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only in isolated points of the gyrokinetic phase-space spanned by the vector (q', 

Then a necessary and sufficient condition that the differential 1-forra dG" = f[dq H can 
be represented in the reduced form 

dG" = f[dq n + f' 2 dq' 2 (25) 
is that /g = /3 (q', £',£>(/>') is a first integral of motion, i.e., there results 

#3(q',£W) = 0. (26) 

PROOF 

Both the necessary and sufficient conditions are trivial. In fact, if /g is a first integral, 
since the Lagrangian 1-form is defined up to an arbitrary gauge it follows 

f' 2 dq' 2 = d(fW 2 ) ~ q'M (q',£W) = 0. (27) 

On the other hand, if up to an arbitrary gauge transformation, the equation f' 2 dq a = holds 
identically in a finite subset [neighborhood] of gyrokinetic phase-space, it follows necessarily 
Eq.(j26]). 

Provided the hypotheses of the lemma hold the following theorem has the flavor of: 
Theorem 1 - Existence of canonical gyrokinetic variables 

In validity of the hypotheses of the Lemma, provided the gyrokinetic transformation 

(rV.jy,^) - (<A<A/i',/ 2 W,0') (28) 

is a -diffeomorphism, it follows that: 

A) it is always possible to identify 

P'l = fi', (29) 
P2 = h'- (30) 

m Eq. H3J); 

B) the transformed variables (q n ,q' 2 ,p[ = fi,P2 = H-,V<t>'-> 0') are canonical gyrokinetic 
variables. 

PROOF 

To prove the theorem one has to realize, first, that the assumptions of the Lemma are 
indeed satisfied by the gyrokinetic Lagrangian defined by EgfTol Then the proof is an 
immediate consequence of the Lemma. 
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A basic consequence of the theorem here pointed out is that the adoption of canonical 
gyrokinetic variables in gyrokinetic theory is only permitted if the gyrokinetic Lagrangian, 
besides ft, has an additional ignorable coordinate, q' 3 and hence it admits necessarily two 
first integrals of motion p'^ and p' 3 . In turn, one can show that this condition implies that 
both the electric and magnetic fields (as well the corresponding EM potentials $, A) must 
be symmetric 2(|. This implies that if the equilibrium magnetic field B is non- symmetric, 
or more generally is locally chaotic (i.e., it does not admit locally a family of nested magnetic 
surfaces), the gyrokinetic transformation / TiPj) - in the sense of approach A - does not exist. 

IV. DISCUSSION AND CONCLUSIONS 

In this paper the conditions of existence of the canonical gyrokinetic variables for a clas- 
sical charged point-particle have been investigated. We have shown that - in the framework 
of an exact gyrokinetic treatment (Approach A) - these variables can only be achieved pro- 
vided the particle gyrokinetic Lagrangian is symmetric. This means, actually, that it must 
have generally two ignorable coordinates (ft and q' 3 ). 

The extension of these results to the asymptotic gyrokinetic treatment (Approach B) 
is non-trivial. In fact, even small perturbations of the EM field can in principle produce 
significant local (and even non-local) stochastic effects. Nevertheless, near an axi-symmetric 
MHD equilibrium, i.e., for magnetic fields which are weakly non-symmetric (and weakly- 
turbulent) - in the sense that they are characterized by suitably small deviations from a 
symmetric equilibrium - one should expect CGKT to hold locally, at least, in an asymptotic 
sense, a result earlier pointed out by White 



However, these conclusions cannot be extended to general situations. As an example, in 
Stellarators magnetic surfaces may only exist locally namely in the neighborhood of nested 
magnetic surfaces only. Therefore it would be highly desirable to be able to construct 
gyrokinetic canonical variables which result independent of the magnetic field geometry and 
apply also to the case of chaotic magnetic fields. An example is given by so-called quasi- 



symmetric 21[ MHD equilibria which arise in Stellarators. These equilibria - which actually 
may be strongly non-symmetric - are expected to be characterized, at most, by a family of 
isolated nested toroidal magnetic surface. Typically, in the intermediate regions between 
these surfaces the magnetic field is chaotic. Another typical situation is that arising in 
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the presence of local MHD/kinetic turbulence, in which EM perturbations may give rise 
to local chaotic behavior of the magnetic field. These results are potentially relevant for 
their implications for theoretical investigations and numerical simulations of magnetized 
plasmas. In fact, the regularity conditions on the EM fields, to be imposed for the validity 
of CGKT, may be locally violated in typical MHD equilibria. For example, a consistent 
kinetic description or a numerical gyrokinetic particle simulation of a magneto-plasma in 
these variables cannot be achieved unless the EM field is weakly non-symmetric in whole the 
domain occupied by the plasma. 

It should be stressed that there is a simple alternative to the description based on canon- 
ical gyrokinetic variables. This is represented by the super-abundant canonical gyrokinetic 



12j , which preserves both the Hamil- 



treatment (super-abundant CGKT) pointed out in Ref. 
tonian character of the equations and - unlike CGKT - is applicable also in the presence of a 
chaotic magnetic field. Basic features of this approach are in fact that: 1) no symmetry (or 
quasi-symmetry) assumption is required for the magnetic field, so that it holds also in the 
case of chaotic magnetic fields ; 2) the Hamilton equations for the canonical pair ( r',p T t) 
are in vector form. Its formulation is summarized by the following constrained variational 
principle, 

Theorem 2 - Superabundant CGKT 

Let x = (r,p r ) be the canonical state of a charged point particle described by the Hamilto- 
nian 



H(r, p r , t) 



1 

2m : 



Pr 



(I 



A 



sc. 



q 



and introduce the diffeomorphism 



(x.' =superabundant canonical gyrokinetic state), where 



1 a * 

Pr' = —A 

EC 



dC 



q 



9 (|r') 



mv'+— A* 

EC 



(31) 



(32) 



(33) 



and the gyrokinetic Hamiltonian (T72j) is represented in the form 



K(^,t) = -p 4>l n , + 7 ^p rl -^-A*(v',u',p^,t) +^V, M ',p^t). (34) 

It follows that: 1) x'(£) is the extremal curve of the functional S'(x') = 
J dt < r • p r i — (p p'^ — K > which satisfies the synchronous variational principle 5S{x!) = 0, 



I 2 



q. 



n 



with the variations #x' = (5r',5p T >, 8(f)', Sp^i) to be taken as linearly independent and the 

function u'{t) to be considered extremal with respect to Sp^i and 5r' , i.e., such that there 
results identically 

5p r , ~ = 0\ (35) 

9p r > 

Sr'-^ = 0; (36) 

2) x'(i) is canonical with respect to the gyrokinetic Hamiltonian 
PROOF 

The proof is straightforward. In particular, by taking the variations with respect to p r > 
and r', the Euler-Lagrange equations for r' and p r > are simply 



d , 1 



-r 



K(*,t) = - p r , -*-A* . (37) 
at op r > ml ec \ 

d Pr , = -l-K(x',t). (38) 



dt dr 

Finally the equations for p', and 0' are manifestly 

& = -±KV,t)=0, (39) 

(40) 

It is immediate to prove that these equations coincide with the equations of motion obtained 
from the Lagrangian fllSH . 
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